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Quantum antiferromagnet at finite temperature: a gauge field approach 

Oleg A. StarykhQ, George Reiter 
Texas Center for Superconductivity, University of Houston, Houston, TX 77204-5932 

Starting from the CP^~^ model description of the thermaUy disordered phase of the D = 2 
quantum antiferromagnet, we examine the interaction of the Schwinger-boson spin- 1/2 mean- field 
excitations with the generated gauge (chirality) fluctuations in the framework of the 1/N expansion. 
This interaction dramatically supresses the one-particle motion, but enhances the staggered static 
susceptibility. This means that actual excitations in the system are represented by the collective 
spin-1 excitations , whereas one-particle excitations disappear from the problem. We also show 
that massive fluctuations of the constraint field are significant for the susceptibility calculations. A 
- - , connection with the problem of a particle in random magnetic field is discussed. 

^ ■ PACS numbers: 75.10; 75.30D 

Some time ago, two self-consistent approaches to the finite temperature properties of the low-dimensional magnets 
were developed: the modified spin wave theory [ |l[ and the Schwinger-boson theory [ ^. Both of them correctly 
describe the generation of the correlaton length by the thermal fluctuations, by complying with the Mermin- Wagner 
theorem [ Due to its conservation of the spin rotational invariance in the thermally (or quantum) disordered 
(paramagnetic) phase the Schwinger-boson theory is a suitable starting point. Moreover, it allows the generalization 
of the Heisenberg Hamiltonian to that of an S'C/(iV)-invariant model [g], for which a powerful 1/N expansion can be 
developed [ Perturbation expansion in a parameter 1/A'^ around the mean- field {N — oo) saddle point provides 
a convinient way of treating the fluctuations in the isotropic system. 
, ■ The applicability of the Schwinger-boson transformation of the spin operator relies completely on the necessary 
constraint at every site: the number of boson degrees of freedom has to be equal to NS, where S is the the value 
and N is the number of boson "colors" [ ^. Due to the constraint only — 1 out of N bosons per site are really 
independent. Usually , this constraint is introduced through the Lagrange field A, defined on every lattice site. In the 
, mean-field approximation on-site constraint is approximated by the averaged one [ ||, and it makes all N bosons free. 

■ As was shown in [ Q , the fluctuations of the constraint field A above the mean- field value and the phases of the 
Hubbard-Stratonovich decoupling field form three components of the C/(l) gauge field in the continuum limit. Also 
it was proved that the continuum limit of the Schwinger-boson large- A^ theory coincides with the same limit of the 
(jpN-i ]-|2odel at distances much larger than the lattice spacing [Q. 

Despite the fact that CP^~^ model at A^ = 2 is just another representation of the 0(3) nonlinear a model, there is 
' no direct correspondence between these models at arbitrary A^ . In principle, identical results should be obtained from 
^ , the exact solutions of both models for the case of physical spins. It is the perturbative expansion around different 

■ N = oo saddle points of the models that makes the predictions of the models very different. Namely, mean-field 
^ \ analysis of the CP^~^ model in the disordered phase produces two branches of elementary excitations with spin- 1/2, 
O • whereas the 0{N) nonlinear a model gives triply degenerate spin-1 excitations in the disordered phase at A^ = cx3 [ |]. 

\ At the same time expressions for the spin-correlation length obtained from these models in the renormalized classical 
^ • region are similar at the mean-field level. 

In this communication we investigate the corrections beyond the mean-field level in the CP^^^ model for the special 
, case of the thermally disordered phase (or renormalized classical phase at T 7^ 0, which corresponds to the Neel-type 
' state at zero temperature) of the GP^~^ model, by exploiting the abovementioned equivalence of it to the continuum 
limit of the Schwinger-boson theory. It is neseccary to stress here that in the analysis to follow we consider the lowest- 
energy excitations with momentum (and energy) much smaller than the inverse correlation length. In that regime 
thermal fluctuations of locally Necl ordered regions destroy completely long range Neel order and antiferromagnet is 
disordered, while inside the regions of the size of correlation length the system behaves as in the broken symmetry 
phase [ ||. The model contains two collective modes: the massless U{1) gauge field (due to the famous mechanism of 
the dynamical generation of the kinetic energy of the gauge field [ ) and the massive fluctuations of the constraint 
field A. A-fluctuations produce only the short-range density-density interaction term, whereas the coupling to the 
gauge field produces the long-range current-current interaction. Due to their masslessness, the gauge fiuctuations 
change significantly the mean-field picture. As we will show here, it is the interaction with the generated gauge 
field that supresses spin-1/2 excitations, while only slightly affecting collective excitations with spin 1. At the same 
time, interaction of the massive A-fluctuations with the mean-field excitations is much weaker, nevertheless it gives a 
susceptibility correction of the same order of magnitude as the interaction with gauge fluctuations. 

Effect of the interaction with generated gauge fleld on the properties of free spin-1/2 spinous was investigated in 
ref.[^ for the case of quantum disordered phase of the SU{N) quantum antiferromagnet, and it was found that due 
to the "hedgehog" -like instanton tunneling events spinous are confined into the spin-1 pairs. Contrast to these we are 
considering the thermally disordered phase without topological defects. 
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Our starting point is the nonlinear a model Euclidian action in the continuum limit 



c r r 1 



25 7o J X^^^^''' + id^nyj,i = x,y (1) 

with the constraint = 1. Vector n describes local staggered magnetization. Putting c — 1, using the CP^ 
representation n = z+crz (z is a complex two-component field), and introducing — —^{z^d^z — {d^z^)z) the 
action becames [ ^ 

S = - f dT I d?v{\{d^~iA^)z\'' + i\{\z\' -I)] ,^l = T,x,y (2) 
9 Jo J 

where Lagrange multiplier A enforces the constraint z^z — 1. 

Generalizing the doublet z to the TV-component vector we arrive at the CP^~^ model action. The saddle point 
equation for the A-field produces the gap m in the spectrum of z-quanta, and finally we get (we have rescaled z by 

V9) 



Z = JdzJdaJ DX'exp- dr J d^r {\{df, - iAi,)z\^ + m^\z\'^ + i\'\z\^} , 
where we have used parametrization [ : 

iX{x, t) = -I- iy{x, r), 



(3) 



with the condition 



J\t J d^rX'ix, 



t) = 0. 



It is important for us that in the thermally disordered phase the gap m is exponentially small as a function of 
temperature [ ^ 

m~Texp (4) 

We have briefly described the derivation of the CP^^^ model from the nonlinear a model for the sake of completeness 
of our discussion. Microscopical derivation of the expression (3) was done in [ ^ , where the coupling of A field to the 
spin-Peierls order parameter was also found , but we will not consider this coupling here. 

Integration over z fields produces the effective action for A and A' fields. In the large- iV limit, the gaussian 
(quadratic) fluctuations around the saddle point <A>=0,<A'>=0 are leading, because the higher orders 
of A (A') will have additional factors [ |^. This quadratic approximation is equivalent to the calculation of the 
polarization operator of z fields. The corresponding diagrams, describing generation of the gauge field energy [|^, are 
shown on Fig.l and in the Coulomb gauge {Aq — 0) an explicit expression for the polarization is 



1 f dp 1 f d^ 

^t.uiq,uJn) = J J^G{p,ei)G{p + q,ei+ujn)i2p + (i)f,{2p + ci)^-26f,^-J2 J 

where G{p, e) is the Green's function of z-bosons in the absense of gauge fields 
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^G{p,ei), (5) 



and ei = 2nlT, ujn = 2'KnT are bosonic Matsubara frequencies. 

Both one-loop integrals in (5) are ultraviolet divergent but the Pauli-Villars regularization makes them finite. 

Because the gap m is much smaller than the temperature T, the lowest excitations (i.e. those inside the gap) 
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correspond to a;„ = in (5), which is equivalent to a static approximation. Making an expansion over ^ in the 
integrand we get 
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Using 

^ ^ e/2 + m2 ^ m^2 ^ e^" - 1 ~ 
we obtain the following expression for the Green's function of the gauge field 



D,M) - - < TA,{q)A,{~q) -^n;J(q) = -^^{5,, - Mfl), (7) 



which describes static random magnetic fields with the correlator 



< h{q)h{-q) ^ (8) 

Physically, the magnetic field h(r) describes the chirality through [ 

]^n-{d^ny. d^n) ^ d^j,A^ - d^A^, ^ h. (9) 

The nature of static approximation is easy to understand: the minimal energy of the z-boson is m, whereas A fields 
are lying inside the gap m and, hence, they can be considered as the static ones in comparison with fast z-quantas. 
Knowing the Green's function, it is easy to calculate the contribution of the gauge fluctuations to the specific heat 

5Cv — m^logT. 

Because of the m? prefactor, this contribution is small. 

The problem of a particle subject to random magnetic field has arisen in the gauge theory of the t — J model [ 
in the vortex lines dynamics of high- Tc materials, and in the study of the transport properties of the normal 
film -superconducting film sandwiches in a magnetic field [ Significant progress has been made recently in 

understanding of the problem, mainly due to the development of nonperturbative approaches [ ^ . Due to the long 
range nature of the vector potential A and two-dimensionality of the problem, the interaction between the 2;-particle 
and the gauge field is strong enough to produce divergences in the leading corrections to the particle's self-energy and 
vertex [||,^. 

As we will show later, z-bosons also have similar divergences. The divergence of the self-energy implies that the 
averaged over the gauge fluctuations Green's function is zero, or the z-particle does not propagate. But it is necessary 
to keep in mind that z-boson describes an = ^ excitation at TV = 2, whereas physical excitations (i.e. those 
corresponding to poles in the susceptibility) are described by the pair of the z-particles , which form an 5 = 1 
excitation [ . 

At the mean field level the susceptibility is given by the bubble diagram (Fig. 2), for an A-component boson system 

X^f{q 0, c = 0) = A-i- ^ = (10) 

X™^((? — Q,uj — Q) diverges at T = (m = 0) indicating the transition into ordered Neel-type state (since n is local 
staggered magnetization, q = Q actually corresponds to antiferromagnetic vector (tt, tt) in the original problem). The 
first 1/A correction to x"^^ is generated by the self-energy insertions to the boson lines and by the vertex correction 
(Fig. 3, solid line denotes z propagator, wiggly - gauge-field propagator). The corresponding analytical expressions 



xy^(g = 0,c. = 0) = ^^ /(^/ (^G3(fc,eOG(fc-g,Q)i?^.(g)(2k-q)^(2k 



xTiq -0,. = 0) = ^J:J^J ^GHk,e,)DM, (11) 
Aj^{q = 0,u; = 0)^^J2 J-^l -^G\k,ei)G'ik-q,ei)D,4q)i2k~ci),i2k-ci). 
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Diagram 36 can be calculated exactly, while expansion over is necessary in other expressions. Straightforward 
calculations give the following answer 



x!'"(, = 0.. = 0) = ^ |i: j^^a j^^Y. (jfiW* /'^'} 

Numerical coefhcients Ci and 5*^ are collected in Table I, as well as symmetry factors (F) of diagrams. 

As one should expect from the gauge invariance, the logarithmically divergent term disappears from final expression 
for the first 1/N correction to susceptibility and, using m as an upper cut-off in the integration over g, we get (subscript 
A denotes the gauge field origin of this correction) 

Xy^(Q = 0,u; = 0) = -^ (13) 

Unexpectedly, the susceptibility correction Xa^ is positive (x = x"^'^ + xY^); which impHes that the gauge 
fluctuations enhance the short range antiferromagnetic order. This result is not as strange as it at first appears. Let 
us recall the phenomenon of the quantum stabilization of the classical order in the Ji — J2 model [ |l3[| : quantum 
fluctuations act to reinforce the classical order beyond the classical domain of stability. 

The finite result for Xa^ means that despite the dramatic supression of the one-particle motion (i.e. the divergency 
of the single-particle self-energy — diagrams 3a and 36), the coherent propagation of a pair of z-particles is not affected 
significantly by the gauge fluctuations. This, in turn, means that the final physical picture is similar to that of the 
0{N) nonlinear a model at iV = cx3 : the elementary excitations in the disordered phase of the antifcrromagnet are 
spin — 1 excitations. 



It is possible to show the origin of the obtained result using an eikonal expansion [ 14 1 . Analogous to [ g] , the 
self-energy divergence corresponds to the vanishing of the z-particle Green's function. In the eikonal approximation, 
the Green's function depends upon the field A through the phase exp{iT) — exp{i (irAp(r)^^), where r(r) is 
a straight line path connecting the initial (r(0) = 0) and final (r(t) = r') points. In the simplest approximation 

r(r) — r'j, hence F = J Iq ^''''t'^tJ-i^)^^'^'^ ^ — ^'ti I (2^)^ ^t^il) small enough q. After averaging over gauge 



field fluctuations we have 



< e'' >= exp--< r- >~ exp--ry^ J < A^{q)A,,{-q) > 



Using eq.(7) we arrive to the logarithmically divergent expression 



< >= exp — -r — — / dOsin 



r N Jo (27r)2g2 

Thus, < e'^ for any nonzero r'. At the same time, in the susceptibility calculations we are dealing with closed 
trajectories [ ^, where by the Stokes theorem / dr A{r)^ is the magnetic flux through the loop formed by trajectory 
r(r) [ 0. According to (8) the magnetic field correlator is constant. Therefore the momentum integration will not 
produce infrared divergency. 

We also tried to calculate the 1/A^ correction to the staggered dynamical susceptibility. In that case, the analogous 
cancelation of the logarithmically divergent terms was found. As before, it is nothing but the consequence of the 
gauge invariance of x- Unfortunately, no damping was found for the frequency less than the gap. Thus, interaction 
with gauge field can not eliminate (at least at 1/iV order) another shortcoming of the mean-field Schwinger-boson 
theory - the absence of the relaxational modes in the system. Such hydrodynamical excitations were found by the 
summation of the ladder diagrams describing the boson-boson interaction beyond the mean-field approximation of 
the Schwinger-boson theory [ |l^. In addition, recent study of the 0{N) quantum nonlinear a model by an 1/A'^ 
expansion [ ^ shows that an interaction with massive fluctuations of the constraint field does produce damping in the 
susceptibility and nonzero density of states inside the gap region. 

Because of all these facts, we understand the problem as following: the main role of the gauge fluctuations is not 
only to eliminate the mean-field excitations from the spectrum, but to conserve the coherent pair propagation. After 
that, the problem is equivalent to the initial nonlinear a model, which naturally describes S — 1 excitations. 

To confirm this point we will calculate the effect of the massive A- fluctuations on x- As in the case of gauge 
fluctuations, we need the polarization operator 11;^ (g, oj) for A-field, which is the dynamical susceptibility near q^AFM = 
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(tt, 7r), (Fig. 2). To be consistent with previous consideration, only static part of Il\{q,uj) has to be kept. The same 
arguments about long-distance behaviour permit us to use an expansion over q^/rri^, and we obtain [|5[| 

0) = E / ^^(P^ + 0(0- = ^^(^. \ (14) 

Correspondingly, the A- field temperature Green's function is -D'^((z) — — pjUx(q o) • Unlike the gauge fluctuations, the 
fluctuations of the constraint field do not produce infrared divergency in the z-particle self-energy. Nevertheless, they 
give correction to x of the same order of magnitute as A-fluctuations. This can be shown by the calculation of the 
diagrams of Fig. 3 , using instead of Z?^^ the A-field Green's function D\ and the corresponding vertex i. By this way 
we obtained (again using m as an upper cut-off in the q-integration) 

i/jv 0.875 

Thus, being relatively unimportant for single-particle properties, fluctuations of the constraint field are as significant 
for the spin susceptibility as those of the gauge field. Summing both corrections together, finally we have 

X™^(g = 0, - 0) + x'^'^iq = 0, = 0) = x"^'(<Z = 0, = 0)(1 + -j^) (16) 

The resulting effect of the fluctuations is to enhance the short range antiferromagnetic order, in agreement with the 
results of [ 1^ . Interestingly, this enhancement arises from the integration over spatial scales which are much longer 
than the correlation length. At the same time, the renormalization of the correlation length and the damping effects 
arise at smaller spatial scale, as was shown in [ ^. 

Equation (16) shows that both types of fluctuations are important in the problem under consideration. To obtain 
the complete 1/N correction to the susceptibility it is necessary to carry out the calculations in the entire energy and 
momentum region, which is much more complicated problem [ ^ and is beyond the scope of our paper. 

In conclusion, we have shown that the spin-1/2 mean-field excitations in the thermally disordered phase of the 
quantum antiferromagnet are suppresed by the interaction with generated gauge massless field. At the same time this 
interaction enhances the susceptibility, which describes the spin-1 collective excitations. Analogous enhancement is 
produced by the massive fluctuations of the constraint field. 



ACKNOWLEDGMENTS 



The authors would like to thank Z.Y.Weng for collaboration on the early stage of this work and helpful discussions, 
and D.N.Sheng and C.S.Ting for useful remarks. We'd like to acknowledge also useful comments of the referee. One of 
us (O.A.S.) is grateful to A.F.Barabanov, A.G.Abanov, V.Brazhnikov, and I.Gruzberg for stimulating discussions. The 
present work is supported by a grant from Robert A.Welch Foundation and the Texas Center for Superconductivity 
at the University of Houston. 



* On leave from the Institute for High Pressure Physics, 142092, Troitsk, Moscow region, Russia. 
^ M.Takahashi, Phys.Rev.Lett. 58, 168 (1987); Phys.Rev.B 36, 3791 (1987); J.Phys.Soc.Jpn. 58, 1524 (1989) 
^D.P.Arovas, A.Auerbach, Phys.Rev.B 38, 316 (1988); A. Auerbach, D.P.Arovas, Phys.Rev.Lett. 61, 617 (1988) 
^ N.D.Mermin, H.Wagner, Phys.Rev.Lett. 17, 1133 (1966) 

"N.Read, S.Sachdev, Phys.Rev.Lett. 62, 1694 (1989); Phys.Rev.B 42, 4568 (1990) 

'^S.Sachdev, J.Ye, Phys.Rev.Lett. 69, 2411 (1992); A.V.Chubukov, S.Sachdev, J.Ye, Yale preprint (1993) 
® A.D'Adda, M.Luscher, P.DiVecchia, Nucl.Phys.B 146, 631 (1978) 

A.M. Poly akov. Gauge Fields and Strings, Harwood, New York (1987) 
® P.A.Lee, N.Nagaosa, Phys.Rev.B 46, 5621 (1992) 
^ L.B.Iofle, V.Kalmeyer, Phys.Rev.B 44, 750 (1991) 
1° L.Brey, H.A.Fertig, Phys.Rev.B 47, 15961 (1993) 

" B.L.Altshuler, L.B.Iofle, Phys.Rev.Lett. 69, 2979 (1992); D.V.Khveshenko, S.V.Meshkov, preprint 1993; J.M.Wheatley, 
T.M.Hong, Phys.Rev.B 43, 6288 (1991); G.Gavazzi, J.M.Wheatley, A.J.Schofield, Phys.Rev.B 47, 15170 (1993) 



5 



A.V.Chubukov, Phys.Rev.B 44, 12318 (1991) 

S.Sachdev, N.Read, Int.J.Mod.Phys.B 5, 219 (1991); F.Mila, D.Poilblanc, C.Bmdcr, Phys.Rev.B 43, 7891 (1991) 
E.S.Fradkin, in "Proceedings of the P.N.Lebedev Physics Institute", vol.29, Consuhants Bureau, New York (1967); 
V.N.Popov Functional Integrals in Quantum Field Theory and Statistical Physics, Reidel (1983) 



TABLE I. Numerical values of C,, Si, and symmetry factors F of diagrams from Fig.3 



Diagram 


C 


S 


F 


a 


1/6 


-1/20 


2 


b 


-1/4 





2 


c 


1/6 


-1/60 


1 
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